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Abstract 

Floods, tides and tsunamis are turbulent, yet conventional 
OQrAodels are based upon depth averaging inviscid irrotational 
^^flow equations. We propose to change the base of such mod- 
^^elling to the Smagorinksi large eddy closure for turbulence 
^^in order to appropriately match the underlying fluid dynam- 
^TiCs. Our approach allows for large changes in fluid depth to 
^^cater for extreme inundations. The key to the analysis un- 
^jd'erlying the approach is to choose surface and bed boundary 

.conditions that accommodate a constant turbulent shear as 
^vga nearly neutral mode. Analysis supported by slow manifold 

theory then constructs a model for the coupled dynamics 
'^«.6f the fluid depth and the mean turbulent lateral velocity. 
l~^The model resolves the internal turbulent shear in the flow 
Q»i}and thus may be used in further work to rationally predict 
(~|erosion and transport in turbulent floods. 

I ^K eywords: turbulent flood, tsunami, Smagorinski clo- 
sure, shallow water equations 
^ ■ 

(^PACS: 47.27.E-, 47.27.nd, 47.10.Fg, 47.11. St, 92.10.Sx 

a^ : 

cnl Introduction 
in 

f~~)This letter models the two dimensional turbulent flow of a 
OOliyer of fluid over flat sloping ground. The fluid of thick- 
^^ness ri(x, t) flows with mean lateral velocity u(x, t). FoUow- 
^ii;ig earlier but more complex modelling using the k-e model 
l^oi turbulence [9j, we approximate the turbulence within the 
rSfl'uid by the nonlinear constitutive relation of the Smagorin- 
j^ski closure [3 HH e.g.]. Then the systematic analysis de- 
'v'eloped in this letter recommends using a nondimensional 
turbulent model described by 
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where Qx and are the nondimensional components of grav- 
ity along and normal to the flat ground, respectively. Fluid 
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is conserved through ^ . The mean momentum equation 
incorporates effects of inertia, Ut, self-advection, uux, bed 
drag, |u|u/ri, gravitational forcing, (g^ — gz^Hx), and an 
enhanced lateral turbulent mixing described by the last 
term. For example, for turbulent flow down an inclined 
flat plate with lateral gravity gx, the nonlinear bed friction 
may balance gravitational forcing whence the above model 
predicts the equilibrium flow to have mean lateral velocity 
u w IS^gxT) . The model (Il])-([2]) also resolves instabilities 
from this equilibrium flow and the emergence and interac- 
tion of solitary or roll waves on the turbulent flow. In a 
personal communication, Howell Peregrine commented that 
environmental shallow water waves were observed to prop- 
agate 1-2% slower than predicted by conventional shallow 
water theory; the above turbulent model successfully pre- 
dicts just such slight slowing because of the slightly reduced 
effect of gravity in the coefficent 0.98 in 

This letter puts the model within the sound sup- 

port of modern dynamical systems theory. Section [31 to em- 
power us to systematically control error, assess domains of 
validity, comprehensively account for further physical effects, 
and resolve internal structures within the turbulent flow. For 
example, Harris et al. [4] modelled particle driven, gravity 
currents using shallow water equations that resolve the dy- 
namics of both the fluid thickness and the mean lateral veloc- 
ity. However, such modelling of turbulent dissipative flows 
from the laminar inviscid foundation of shallow water equa- 
tions appears a contradiction that demands resolution using 
the sort of approach introduced here. Future developments 
using this approach may straightforwardly incorporate com- 
plex physical effects such as the three dimensional flow over 
varying terrain similar to previous modelling of laminar vis- 
cous flow j8j. 

2 Differential equations of the 
Smagorinski large eddy closure 

Let the incompressible fluid have thickness r|(x, t) above the 
ground, constant density p, and flow with some turbulent 
mean velocity field u = (u,w) = (ui.ua) and turbulent 
mean pressure field p. In this letter we restrict attention to 
two dimensional turbulent flow. 

Smagorinski large eddy closure Define the turbulent 
mean strain-rate tensor [TS^, e.g.] 
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where X] = x and X3 = z are distances along and normal 
to the mean solid substrate, respectively. Then take the 
turbulent mean stress tensor for the fluid to be (Jij — — pSij + 
IpwEij for some effective turbulent eddy viscosity y. The 
Smagorinksi model [T^i e.g.] for very large Reynolds number 
flows then sets the turbulent eddy viscosity "v to be linear 
in the magnitude £ of the second invariant of the strain-rate 
tensor: 

■V := CtTi^e where := ^ . (4) 

i-J 

Comparison with established channel flow experiments 
e.g.]. Section [O] indicates the constant of proportionality 
Ct ~ 0.02 for the turbulent environmental flows we consider. 

The lateral turbulent mixing in the model strength 
0.26riu, is an order of magnitude larger than direct mixing in 
the underlying Smagorinski model, "Vt « 0.02ri^-|^ , because 
of the 'shear dispersion' enhancement to lateral mixing via 
the interaction between lateral shear and vertical mixing. 
Depth averaging misses this important physical interaction. 
This enhancement in ^ qualitatively agrees with Bijvelds 
et al. jl] who increased the lateral eddy viscosity in order for 
their numerical model to match laboratory experiments. 

Partial diflferential equations Make variables nondi- 
mensional with respect to some velocity scale, a typical 
fluid thickness, and the fluid density. Then the nondimen- 
sional pdes for the incompressible, two dimensional, turbu- 
lent mean, fluid flow are firstly the continuity equation 

and secondly the momentum equation 

^+u- Vu = -Vp-FV-T-Fg, (6) 
at 

where g — (g^, Qz) is the nondimcnsional forcing of gravity, 
and T is the nondimcnsional, turbulent mean, deviatoric, 
stress tensor0 

r,.2.m,^c.,H{^^^. (7) 

Boundary conditions Compatible with the Smagorinski 
closure, wc formulate boundary conditions on the ground in 
terms of the turbulent mean velocity and the fluid depth. 
On the mean bed allow slip to account for a relatively thin 
turbulent boundary layer, 

w = and u = Cu'n^ on z = , (8) 
oz 

where comparison with experiments on channel flow. Sec- 
tion |4Tl indicates the coefficient « 1.85; in applications 
this coefficient would vary to parametrise different ground 
roughness. We envisage that z = denotes the local mean 
level of the ground in analogy to u and p denoting the tur- 
bulent mean velocity and pressure fields. 

-"^Some [51 e.g.] like to use the negative of this stress tensor, but then 
"+V • t" in JSJ has to be "—V ■ t" . 



The kinematic condition on the turbulent mean free sur- 
face is 

— +u—^v onz = ri, (9) 
at ox 

where similarly r\ denotes the turbulent mean location of the 
free surface. 

The turbulent mean stress normal to the free surface 
comes from constant environmental pressure, taken to be 
zero, that is, 

-P+ ,J 2 (^33 -2ti^Ti3 -m^Tii) ^ onz = ri. (10) 

Lastly, there must be no turbulent mean, tangential, stress 
at the free surface: 

(1 -ri^)Ti3 -hTlx(T33 -Til) =0 0nZ = Tl. (11) 

3 Centre manifold theory supports 
the modeUing 

This section describes the new approach to placing mod- 
els such as (II])-© on a sound theoretical base. We aim the 
mathematics to represent the dominant physics: as indicated 
by Figures 14-15 by Janosi [5], turbulence does mix across 
the fluid layer. Thus we expect to be able to model the dy- 
namics in terms of depth averaged quantities. But instead of 
crudely depth averaging the equations, we use centre man- 
ifold theory to resolve turbulent dynamics across the fluid 
layer and hence provide a sound macroscale closure for the 
relatively slow, long term, dynamics of interest to environ- 
mental modellers. 

To apply the theory we artiflcially modify the zero tangen- 
tial stress free surface condition (fTT|) to have an artiflcial forc- 
ing proportional to the local velocity, a forcing parametrised 
by Y and which we later remove by evaluating at parameter 
Y = 1 in order to recover the physical turbulent equations: 
on z = T) replace (fTTj) with 

[(1-Tl"jTi3+T1x(T33-Tii)] ^ j} ^ j\ u\ (12) 

v2( I -I- Cur 

Evaluated at y = 1 this artificial right-hand side becomes 
zero so the boundary condition reduces to the phys- 
ical boundary condition of zero, turbulent mean, tangen- 
tial stress (fTTj) . However, when the parameter y = and 
the lateral gravity and lateral derivatives are negligible, 
Qx = 9x = , a neutral mode of the dynamics is the lat- 
eral shear flow u (x Cu + z/ti . This neutral lateral shear 
mode arises because in pure shear flow the stress component 
Ti3 = ydu/dz and hence the artiflcial free surface condi- 
tion reduces to -vdu/dz = "vu/ri on z = r| . 

Conservation of fluid provides a second neutral mode in 
the dynamics. That is, when y = gx = 9x = , then a two 
parameter family of equilibria exists corresponding to some 
uniform lateral shear, turbulent mean, flow, u cx Cu+z/t) , on 
a fluid of any constant thickness r|. For large enough lateral 
length scales, these equilibria occur independently at each 
location x [T31II11 e.g.] and hence the space of equilibria are 
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in effect parametrised by u(x) and ri(x). Provided we can 
treat lateral derivatives 9x as a modifying influence, that is 
provided solutions vary slowly enough in x, centre manifold 
theorems [31 [71 e.g.] assure us of three vitally important 
properties: 

1. this subspace of equilibria is perturbed to a slow man- 
ifold, on which the evolution is slow, that exists for a 
finite range of gradients 9x, and parameters y and gx, 
and which may be parametrised by the average, turbu- 
lent mean, lateral velcity u(x, t) and the local thickness 
of the fluid r|(x, t); 

2. the slow manifold attracts solutions from all nearby ini- 
tial conditions exponentially quickly on a cross-depth 
mixing time; and that 

3. a formal power series in the parameters y, Qx and gra- 
dients 9x approximates the slow manifold to the same 
order of error as the order of the residuals of the gov- 
erning differential equations. 

That is, the theorems support the existence, relevance and 
construction of slow manifold models such as id])-©. 

An alternative and powerful view of these theorems is that 
they follow from a nonlinear, normal form, coordinate trans- 
form that decouples the slow and fast modes in the fluid 
dynamics [T31 e.g.]. That is, the models we discuss are es- 
sentially simply a reparametrisation of the state space, but 
only that part of the state space where the relatively slow 
dynamics occurs. 

4 Low order models of the dynamics 

The detailed and complicated algebra deriving the model 
is of little interest to users of the model. Computer alge- 
bra readily constructs our slow manifold models [T7]: those 
interested could check the code and verify that the algo- 
rithm asymptotically solves the governing differential equa- 
tions and boundary conditions as specified pp. 9-16]. 
The asymptotic construction of the slow manifold is valid 
for small lateral derivatives, small lateral forcing and small 
perturbation of the free surface condition [Property [3]. This 
section focusses on the resulting model and its interpreta- 
tion. 

4.1 Channel flow calibrates the model 

Open channel flow is well studied. Experiments give em- 
pirical structure of the turbulence in the interior of the flow 
which we use to calibrate the two constants in our Smagorin- 
ski model of turbulent floods. 

The cross fluid structure of the velocity, pressure and 
strain rate forms the slow manifold, here with no lateral 
variations, 9x = , as we compare with uniform flow. In 
terms of the scaled normal coordinate C — z/r\ , computer 
algebra [T7] deduces the lateral shear velocity profile 
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Figure 1: approximate vertical profiles at open channel flow 
equilibrium for which g^r]/u^ = 0.0031 . 
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the magnitude of the rate of strain tensor is 
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and flows with hydrostatic pressure, p — gzTl(1 — C) . These 
fields approximate out-of-equilibrium turbulent flow. The 
corresponding evolution of the mean lateral velocity is 



du 
dt 
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Obtain a physical model by setting the artificial parameter 
Y = 1 . Lateral gravitational forcing, the second term, accel- 
erates the flow until balanced by dissipation, the flrst term, 
to reach the well known equilibrium for open channel flow. 

Observations calibrate the Smagorinski coefficient Ct and 
bed 'slip' coefficient Cu- Many experiments confirm the 
equilibrium flow velocity is ii = 1 8u, for 'friction velocity' 
'U-* '■— \/9xTl • Summarising many experiments, Nezu gave 
empirical formulae for the structure of turbulent energy k 
and dissipation e [H] (17) and (21)]; we integrate the eddy 
viscosity ■y{z) — C^k^/e obtained from these formulae over 
the depth gives the empirical, depth averaged, eddy viscosity 
'V — 0.0796 riu* . To match the equilibrium flow velocity and 
the average eddy viscosity we find Ct = 1 /50 and Cu — 11/6 
which we use hereaftero Figure [T] plots the corresponding 
vertical profiles of the rate-of-strain i and lateral shear ve- 
locity u at the equilibrium of open channel flow: the lateral 
velocity shows a small shear corresponding to the turbu- 
lent mixing across the fluid with dissipation at the bed; the 
rate-of-strain increases towards the bed (imtil the bottom 
boundary layer which is not resolved in Figure [ij. Our reso- 
lution of these internal fields could contribute to prediction 
of erosion and deposition in further modelling. 



^The equations determining these constants appear to be a reason- 
ably well conditioned. 
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4.2 Turbulent spatio-temporal dynamics 



The previous subsection briefly explored the equihbrium of 
uniform open channel flow. Now we address modelling flood 
flow with mean turbulent properties varying in space and 
time. 

Further computer algebra [17] generates modiflcations to 
cross fluid structures such as p^ - lfTS|) to account for rel- 
atively slow variations in the lateral direction x, small but 
non-zero dx- This letter is no place to record the consider- 
able algebraic detail. However, the result is that a model for 
the evolution of mean lateral velocity ii(x, t) is 



f = - 1 .045 - 0.0031 1 + 0.985 f g. - g.^ 

ot ox 11 V ^'^ 



+ 0.058T1 ( ^ 



.0.259n|uS + 0.522 lul^l^ 
ox-^ ox ox 



- 0.015^u2 - O.OOZ^u^ - 0.007u2 1^ 
Ti r| ox^ 

+ 0(Y3/2 + g3/2 + g3 + 93)_ (Ig) 

This would be solved with the associated equation ((T|) rep- 
resenting conservation of fluid. Obtain the model flj-lll]) 
discussed in the Introduction of this article by four steps: 
evaluate (fT6|) at the physically relevant y = 1 to remove 
the artiflce in the surface boundary condition (jl2p : approxi- 
mate coefficients to two decimal places; neglect the numeri- 
cally small effects recorded in the third line of (fTH)) : and, to 
small numerical errors, combine the three terms on the sec- 
ond line of (I16p into the one nonlinear lateral mixing term 
on the second line of ^ . Importantly, although the models 
(HI) and are expressed in terms of depth averaged lateral 
velocity, remember that they are derived, not by depth av- 
eraging, but instead by systematically accounting for inter- 
actions between vertical profiles of velocity/stress and bed 
drag and lateral space variations. Centre manifold theory 
provides the framework to account for these physical pro- 
cesses and so deduce these physical models. In particular, 
this approach provides a sound lateral regularising dissipa- 
tion in the second lines of ^ and (fTH|) . 

Computer algebra [17l §5] constructs terms in the formal 
power series solutions in the notionally small parameters y, 
gx, Qz and dx- We record one truncation of the formal 
power series in (|16p . Other truncations of the multivariate 
power series generate other valid approximations of varying 
accuracy. With the support of centre manifold theory, re- 
searchers may choose an approximate model that suits the 
parameter regime of their application. 

4.3 Turbulent roll waves arise 

Roll waves spontaneously arising on fluid flowing down a 
slope have been noted and analysed for a century [101 [TJ 
e.g.]. Figure [2] plots the proflle of one period of such a turbu- 
lent roll wave on a slope as predicted by the model (IT|)-([2]). 
The occurrence of such roll waves is of concern for hydraulic 
engineers: the model (IT])-([2]) predicts turbulent roll waves 
arise on slopes greater than 1.38%. 




Figure 2: turbulent roll wave on a 4% slope on a fluid 
of nondimensional mean depth 1 and with nondimensional 
mean lateral flow velocity 3.55. The fluid depth (blue) varies 
by nearly 50% and is closely in phase with the velocity vari- 
ations (green). 
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Figure 3: critical slope gx for the appearance of turbulent 
roll waves of wavenumber k as predicted by model 



4 



Standard stability analysis on a fluid of nondiniensional 
depth 1 and the corresponding nondiniensional equilibrium 
flow velocity 1 7.8^g7 seeks small perturbations e exp(ik,x + 
At) . For model (P)-© the characteristic equation for such 
perturbations is 

- (0.1102 + 36.5ik + 4.6231c^ l^g^A 
+ [2.94gxik+ (0.98 -331.9g^)k2+ 82.1 9g^ik^] =0. 

Consequently, the real part of the growth rate A crosses zero 
for critical slopes 



0.00296 + 0.248k2 + 5.2k4 
0.2146 -k2 



(17) 



as plotted in Figure [3] (note that the nondiniensional gx ef- 
fectively denotes the slope of the bed). Roll wave instability 
arises first for long waves on slopes steeper than 1.38%. On 
a slope of 4%, as simulated in Figured! roll waves are unsta- 
ble for nondimensional wavenumber |k| < 0.12, that is, for 
wavelengths longer than 52t). These predictions are experi- 
mentally verifiable. 

5 Conclusion 

Following similar modelling for viscous thin fluid films [181 
[TB] , the innovation of modifying the free surface condition 
to (fT^ with the Smagorinski closure places the modelling of 
turbulent fluid flows upon the powerful and sound basis of 
centre manifold theory [SI [7l e.g.]. Such a modern dynami- 
cal system foundation empowers us to systematically derive 
the models ([1]), ^ and for turbulent tides, floods and 
tsunamis. Earlier modelling based upon the k-e turbulence 
model led to a description of floods in terms of the depth 
averaged velocity, turbulent energy and turbulent dissipa- 
tion [9] . The extra complexity of these extra field variables 
makes the earlier model more problematical to use in fore- 
casting environmental flows. The simpler model here should 
prove more useful. 

The derivation resolves vertical profiles within the flow, 
and the physical processes and interactions between these 
profiles and lateral variations. Thus these models, via 
(fO|) and (fH|) for example, can be used to rationally pre- 
dict bed strain and hence erosion and transport in detail 
for out of equilibrium flows. Further refinement of the com- 
puter algebra derivation 17J will empower sound modelling 
of three dimensional turbulent flows, via the Smagorinski 
closure, on varying topography as already done for viscous 
fluid films [18]. 
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